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Abstract 

^: 

(~| I We study the potential of radiative decays of the T15 and of the (j) niesons to search for a light 
pseudoscalar Higgs boson, proposed as a possible interpretation of S"*" pfi'^fi~ events observed 
I . by the HyperCP collaboration at Fermilab. We conclude that the detection of this signal should 
^ I certainly be possible with the current CLEG T15 data, and is within the reach of KLOE in at least 

' part of the range of couplings suggested by the HyperCP findings. 
I> ■ 

■ 1 Introduction 

o ; 

. ^ ^ ^ 

. The HyperCP collaboration at Fermilab has recently reported [1] the observation of 3 S+ PfJ- IJ- 
^ I decays, with the dimuon invariant mass M^^ = 214.3 it 0.5 MeV, consistent with the decay of a 
narrow neutral intermediate state P^. Among the several candidate new particles put forward to 
explain this observation [2]- [7], the possibility [7] that is a light pseudoscalar Higgs (a^) in the 
5^ \ next-to-minimal supersymmetric Standard Model (NMSSM) interestingly fits within theoretical 
frameworks recently explored in the attempt to loosen the very tight constraints set in the MSSM 
by the non-observation of the Higgs at LEP2. 

Comprehensive studies of the implications of a light NMSSM Higgs for current and future 
measurements have been performed. Excellent examples are given, for example, in [HI [9] and in 
a comprehensive review of non-standard Higgs scenarios [10]. Here we limit ourselves to explore 
one important direct implication of the proposal made in [7]. Namely the possibility that be 
produced in quarkonium radiative decays V ja^ [TI], where, in our analysis, V = Tis or 
(p{1020). The most recent and precise limit on T15 — > 70*^ was given by CLEG [12j, where however 
a" is required to be stable and non-interacting. While several studies with visible Higgs final states 
have been performed since the early works of [13], we are not aware of limits valid in the case 
of prompt — > decays in the range of branching ratios B(Ti5 l^P) relevant to the 

model of [3- Likewise, we are not aware of experimental studies on the potential of (j) decays to 
detect or constrain such scenarios. The excellent performance of the KLGE detector at DA<I>NE 



X 
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(see e.g. [H]), and the prospects for a continued operation at much higher luminositie^, make the 
consideration of (p decays extremely interesting. 

Radiative quarkonium decays have been studied already in the context of other hypothetical 
interpretations of the HyperCP findings, namely the identification of with a sgoldstino [6j. In 
that case decay rates are typically beyond the present experimental sensitivities; in the case of a 
pseudoscalar Higgs, as we shall show, decay rates consistent with the HyperCP data are well within 
the reach of current experiments. 

The virtue of our analysis is that it is only based on the assumptions on the coupling of 
to down- type quarks and to leptons, and it is independent of the details of the pattern of Higgs 
expectation values and mixings which need to be fixed to study other possible implications of this 
class of models [HI [9]. Furthermore, the quarkonium decay does not depend on the size of FCNC 
couplings of a^, which enter in the analysis of the S+ pfi~^fj.~ decay. Finally, the prediction for 
(j) decays only depends on the coupling to the s quark, which is what is probed in S+ — > pfi^fi~, 
and therefore evades the potential loophole of the T15 study, namely that in more complex Higgs 
sectors the coupling to the b quark could be independent of the coupling to the s. Therefore, 
independently of the confirmation or invalidation of the hypothesis of [7], our note should be taken 
as an encouragement to revive the exploration of radiative quarkonium decays as a powerful tool to 
detect or constrain non-standard Higgs models, extending the considerations presented in ref. [8]. 
Other motivations to continue the experimental study of quarkonium radiative decays as a probe 
for new physics are presented in [15j . 

We start by briefiy summarizing the key couplings of to the fermion sector. 



ff = - {gu rriu wy^u + ga ma d-f^d + ga me i'y^i) — , (1) 



V 



where v ~ 246 GeV, and the couplings gu^d depend on the details of the SUSY and EW symmetry 
breaking pattern. The analysis the HyperCP data in [7j suggests gd = 0{1), compatible with the 
constraint \gd\ ^ 1.2 from the muon g — 2[2\. A rough lower limit on gd can be obtained by requiring 
the Higgs lifetime to be short enough to guarantee the reconstruction of a common vertex for the 
three tracks relative to the muons and the proton in the HyperCP data. Using 

r(aO^^V-) = |^ hlo-^mlf' (2) 

and using a Lorentz 7 factor of about 230 (the typical momentum of the dimuon pair in the 
HyperCP data is 50 GeV [20]), the average decay length of the in the laboratory frame turns 
out to be ~ 0.02 cm/g'^. Requiring this not to exceed a resolution in z of the order of 60 cm [2D], 
would lead to gd ^ 0.02. 

The above couplings induce quarkonium decays with a rate given by: 

where mg is the heavy quark mass and F ~ 0.5 includes the effect of QCD radiative corrections. 
For our numerical analysis, in the case of the T15 decays we'll use = Mis/ 2. In the case of 



^See ^http: / /mra .Inf . infn. it /lnf admin /direzione/roadmap/roadmap.htmi] 
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(j), 213.8 MeV< M^^ < 214.8 MeV 
S/VB for 2.5 fb^^ 


0.5 pb 
1-8 9j 


214.2 MeV< M^^ < 214.4 MeV 
S/VB for 2.5 fb-i 


0.1 pb 
3-9 gj 


Ti5, 2m^ < M^^ < 220 MeV 
for 1 fb-i 


0.05 pb 
150 5^ 



Table 1: Cross section for the QED process e+e ^ fi^fj, 7, and significance of the a" 
signal. 

the (j), the small mass leads to a rather large uncertainty on the choice of the quark mass and of 
the size of QCD radiative corrections. As in the case of the T15, we shall adopt the convention 
rus = and discuss in appendix [B] the motivations and limitations of this choice. With these 

assumptions, and using B(Tis e+e") = 2.4 x 10~^ and B^cp — > e~^e~) = 3 x 10~^, we obtain: 

B{Tis^ia^) = l.9xlO-^Fgj , B{ct> ^ ja"^) = 2.8 x IQ-^F gj (4) 

The rate of Tig produced at the B factories depends on the details of the beam energy resolution, 
which is typically larger than the intrinsic T15 width. Using the CLEO data discussed in [16] as 
a benchmark, we extract a production rate of 21 x 10^/1. 13fb^^ ~ 1.8 x lO'^fb. A similar rate is 
expected for BaBar and Belle, due to comparable beam energy profiles. The DA<I>NE beams allow 
instead to operate at the peak of the (p resonance, with a rate given by: 

127r 

a(e+e- <p) = — B,+,- Aisr ~ 2.8 x 10^ fb , (5) 

where Aisr ~ 2/3 is the QED ISR correction factor. This leads to the following estimates for the 
number of F — > 70*^ — > "ffi'^fi^ decays per fb~^ (assuming F = 1/2): 

iV(Tis ^ 7«° ^ = 1700 5d (6) 

N{<p^ja^ ^-ffi+fi-) = 40 gj. (7) 

Considering that the current KLOE dataset consists of about 2.5fb~^, and that a factor of > 10 
increase is planned for the future, we conclude that both the B factories and DA<I>NE are potentially 
sensitive to the range of couplings required for the Higgs interpretation of the HyperCP events. 

The actual observability of this signal depends however on the background rates. While the 
signature is rather sharp (a monochromatic photon, produced at large angle with energy close to 
mv/2 and recoiling against a low- mass dimuon pair), a potentially large background can arise from 
the continuum QED production of e^e^ 7/2^^/2" events. The actual size of the background will 
depend on the dimuon mass resolution, where the signal has a sharp peak with an intrinsic width 
below the MeV. We point out that, for a muon pair produced near threshold, like in the case of 
the HyperCP signal, the invariant mass measurement should have an improved sensitivity. The 
invariant mass (as a function of the angular separation and the momentum of the muons) is in 
fact near a minimum, and therefore has a reduced sensitivity to the measured parameters and their 
errors, as discussed in an AppendixO For our numerical study we assume dimuon mass resolutions 
equal to ±5 MeV for the Tis, and ±0.5 MeV or ±0.1 MeV for the (j). 
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The resulting background rates are given in table [H together with an estimate of the statistical 
significance of a possible signal. Here we required the photon to have energy larger than 1 GeV 
(100 MeV) for the T15 ((/>), and to be emitted at an angle with respect to the beam such that 
I cos^l < 0.7. Since the photon distribution in the signal is proportional to 1 + cos^ 9, this cut has 
a signal efficiency of about O.GlH. 

The separation of the signal from the background appears rather clear at the Ti^- with lfb~^, 
even if the coupling is significantly smaller than its preferred 0{1) value. The case for KLOE 
with its current statistics is at the borderline, even though there is a margin of uncertainty due to 
the large size of the QCD corrections. A mass resolution of 0.1 MeV appears required to yield a 
significant signal. 

Higher statistics, particularly at KLOE, would allow to strenghten the constraints on similar 
explorations of exotic BSM scenarios. 

Acknowledgments We thank D. Asner, S. Bertolucci, F. Bossi, P. Dauncey, G. Isidori and H. 
Park for useful discussions and contributions to this study. 



A Background estimates 

The exact matrix element for the process e^(pi)e+(p2) 7(^)m (91)^^(92) can be easily calcu- 
lated. We show here the result in the case of initial-state emission only, which is the relevant one 
for our signal: 

((piQ)' + (P2Q)')] (8) 

where qi + q2 = Q- The final-state kinematics of the signal, two low- mass muons recoiling against 
the photon, allows us to approximate the cross section as: 

da(7/i+^-) _ dajjj) dQ^ a 
dcosO^ ~ dcosO^ Q2 2tt ^^'"^^ ' 

where 



and 



d cos 0.y ut 32-irs ' 



2 



P{z) = [z' + {l-zf + 2^] (11) 



is the Altarelli-Parisi splitting function, for 



1 / 4m2 1 / 4ml 



^We verified that S/^/B is rather constant, and maximal, when the cut on cos6' is varied in the region 0.6-0.8. 
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After integrating over z, over the dimuon mass in the range 2m^ + < Q < + A2 (with 
Ai^2 << ^fi), and over the photon angular distribution in the range — c < cos 9 < c, we obtain: 

While the results of table 1 were obtained using the exact formula, eq. (fT3]l is an excellent approx- 
imation to the exact result, which can be used to easily explore the background size under cuts 
different from those used in table 1. 



B QCD corrections 



The QCD corrections to the F — > P + 7 decay of quarkonium are given by [T9l 



TiV ^ e+e-) ~ 



(1 - mp/niy) 1 - asCf/n 



ap 



pmax 
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(14) 



where we have assumed for simplicity gd = 1- The —1/4 term accounts for the QCD corrections 
to the e~^e~ decay rate. This result is obtained by renormalizing the quark mass in the on-shell 
scheme. Thus, the quark mass appearing in Wilczek formula is the pole mass; it can be given in 
terms of the MS mass evaluated at the scale of the quark mass itself as |17j 



so that eq. (fT^ becomes 



T{V e+e-) 



vra 



rrin 



-(1 



q{rnq) [1 + C/as/vr] 



m\lra\r) 



a 



TT 



ap 



^max 



(15) 



(16) 



where ap is an increasing function of its argument, with a maximum value ap(l) ~ 6.62. The 
strange MS mass is determined to be around 0.1 GeV at 2 GeV [TT]. Its value at a scale equal to 
the mass itself can be obtained using the renormalization group equation relation 



4/9 



(17) 



and since as{fJp) diverges for fj, < 0.5 GeV, one would always determine fhsifhs) > 0.5 GeV. 
Using 4-loop running for the strange mass [18], as{Mz) = 0.119 and ms(2GeV) = O.lGeV, we 
find ifis{rhs) = 0.6GeV. The next question is what value of Us one should use in eq. (jl6p . The 
coefficient of Ug is ~ 1.85, and the scale of Ug should be taken around the annihilation energy. 
We find as(lGeV) w 0.5, so that at this scale the radiative corrections has the same size of the 
Born term. Assuming that the correction exponentiates, we get exp(— 0.5 x 1.85) ~ 0.4. Assuming 
instead a rational form 1/(1 + 0.5 x 1.85) ~ 0.5. We thus conclude that setting the strange mass 
to half of the (j) mass, and using a correction factor of order of a half yields a reasonable estimate. 
We should not forget, however, that uncertainties are very large, so that a rate a few times larger 
or smaller than our estimate is not unconceivable. 
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C kinematics near threshold 

The invariant mass of the two muons, as a function of their momenta pi, p2 and thair angular 
separation 9 is 

M^fi = + 2pip2{l -cos 9) +2[EiE2-piP2- ml] 



4ml + 2P1P2(1 - cos 9) + 2ml ^ J ^1 2 ■ (l^) 



Near threshold, the second and third terms in eq. (jlSp become small, so that in this limit 

M„,.2™„ + j^«^ + !^^, (1.) 

where p is the common value of pi and p2, E is their energy, and Ap = pi — P2 is their small 
momentum difference. Experimental errors S9 and 5Ap on the measurement of the angle and the 
momentum affect the invariant mass measurement by 

SM^^^^em + '^pA,. (20) 

Thus the error is reduced by the small factors 9 and Ap. 
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